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ABSTRACT 

As neutron star - black hole binaries are expected to be discovered through fu- 
ture pulsar surveys using upcoming facilities, it is necessary to understand various 
observable properties of such systems. In the present work, we study the advance of 
the periastron of such binaries under the post-Newtonian formalism over a wide range 
of parameters. We find that the first and second order post-Newtonian effects and the 
leading order spin-orbit coupling effects arc significant for such binaries but higher 
order effects can be neglected. 



Key words: 

relativity 



stars: neutron — pulsars: general — binaries: close — gravitation 



1 INTRODUCTION 

Relativistic binary pulsars are very useful tools as they 
lead to the test of general relativity, accurate determi- 
nation of the masses of the members of th e binary, un- 
derst anding of the stellar evolution, etc (|Stairs! I2003L 
2004). So far, such systems include only neutron star 
- white dwarf and neutron star - neutron star sys- 
tems. No neutron star - stellar mass black hole (NSBH) 
system is know n yet, although a number of theoreti- 
cal studies exist JPfahl. Podsiadlowski. fc Rappaportl [20051 : 
iKiel fc Hurlevll2009l : lFaucher-Giguere fc LoeblboilT and ref- 
erences therein), discussing the expected number and prop- 
erties of NSBH binaries formed either through normal evo- 
lutionary channel or through exchange interactions in dense 
stellar environments. NSBH systems are expected to be dis- 
covered through future radio pulsar surveys using upcoming 
facilities like SKA. Moreover, gravitational waves emitted by 
these systems will be detected soon by the ground based de- 
tectors like LIGO, VIRGO, TAMA300, and GEO600. For all 
these reasons, theoretical efforts to understand the proper- 
ties of NSBH systems in better details is ongoing. These sys- 
tems are usually modelled under the Post-Newtonian (PN) 
formalism. In the present work, we explore which terms in 
this formalism are more important and deserve more atten- 
tion. We confine ourselves with only one parameter, namely 
the advance of the periastron of the orbit of the neutron 
star, which is a measurable quantity when the neutron star 
is a radio pulsar. 
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The rest of the paper is organized as follows, in Sec- 
tion [2] we present our study, then in Section [3] we discuss 
some other relativistic effects in such systems, and finally 
our conclusions are in Section [3] 



2 ANALYSIS 

PN approximation is effectively an expansion of Einstein's 
theory in powers of a small parameter e which becomes 
(v/c) 2 ~ GM/(c 2 R) in case of a compact binary where v, 
M, R, G, and c are the characteristic orbital velocity, the 
total mass, the typical orbital separation, the gravitational 
constant, and the speed of light. In addition to conventional 
PN parameters, spins of the members of the binary also con- 
tribute to the dynamics of the system and the total Hamil- 
tonian becomes 



H 



II i £jLO | rrLO 
iipN + -"SO + tlSi 



«2 



i rrLO , rrLO , rrNLO , rrNLO 
+ ii, 9 2 + -«.<? 2 +-«SO + H , 



rrNLO , rrNLO , 

-O o2 -\~ -tl e2 + ■ 



S 2 TiJSO -r^S]^ 

higher order spin terms. (1) 



where PN stands for the Post-Newtonian expansion, nth 
order term involving (v/c) 2n . SO is for the spin-orbit cou- 
pling, S1S2, S 2 and S 2 are for the spin-spin interactions. 
LO means the leading order, NLO stands for the next lead- 
ing order, NNLO means the next-to-next leading order and 
so on. The Hamiltonia n is known upto 3.5PN and NNLO- 
SO and NNLO-S1 & (TTessmer Hartung, fc Schlferl l20ld: 



Konigsdorffcr, Favc, & Schafcr 2003; lHartung fc Steinhofli 
2011al lbl: ISteinhoff. Hergt. fc Schaferl |200S| , and references 
therein), but the observable parameters, e.g., the rate 
of the periastron advance of the binary are known at 
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most upto 3PN and LO-SO IPamour fc Schaferl 1 19881 : 



Kb'nigsdorffcr & Gopakumar 2005). 

iKonigsdorff er fe Gopakumar! |2005i . hereafter KG05) 
parametrized the orbital motion of a compact binary up 
to 3PN and LO-SO, with a parameter k which becomes 
the measure of the periastron advance of the binary, i.e. 
k — k = — (u)}, for the non-spinning part. Here Co is the 
rate of the periastron advance and n is the mean orbital 
motion which can be expressed in terms of the reduced en- 
ergy (E), the reduced orbital angular momentum (L), and 
r\ — mim,2/(mi + 1712) 2 , mi and m.2 being the masses of 
the members of the binary (see Eqn. 5.6c of KG05). But 
fcgo 7^ ^so m their formalism. Following KG05, k up to 
3PN can be written as: 



(2) 



2 L 2> 



and 



fepN = &3PN = 



3 -2E 
c 2 L 2 4c 2 



3 (-2E) 2 



(-5 + 2 V + (35 - W V )x) , (3) 



120 - I2O77 + 96J7 2 + 



c 2 L 2 384c 4 

- 1237 

(36960 - 40000?? + 615tt 2 77 + 1680?7 2 )3; 2 ] 



(-10080 + 13952?? - 123rr 2 77 - 1440r; 2 )a; + 



(4) 



where x = l/(—2EL ) and fcpN = &ipn + fepN + &3PN 
(neglecting higher order PN terms). In Newtonian dynam- 
ics, x = 1/(1 — (?) where e is the eccentricity of the 
orbit givin g x ^ 1. The 3PN expressions for E and L 
are given bv lMemmesheimer. Gopakumar. fc Schaferl (j2004T ) 
in terms of different orbital elements, and here it is pos- 
sible for x to be smaller than 1, but it can have only 
such values that the square of the eccentricities are ei- 
ther zero or positive numbers (Eqns . 20b and 20d of 
iMemmesheimer. Gopakumar. fc Schaferl |2004l )). We have 
checked for a wide range of values of (— 2_E)/(4c 2 ) that the 
square of the "radial eccentricity" (e 2 ) remains greater than 
or equal to zero for x > 0.59. So we use only the values of x 
greater than or equal to 0.59. 

Eqn. ((3]) implies that fopN is zero at x = (5 — 2-q) / (35 — 
IO77). For realistic ranges of 77 for NSBH systems, this hap- 
pens around x = 0.14. Here by 'realistic ranges of rf we mean 
realistic ranges of masses for neutron stars (m ns = 1 — 2 Mq) 
and black holes (mth = 5 — 15 Moj]]. In this range, we 
choose three representative values of 77, e.g., 77 = 0.107725 
for mi = m ns = 1.4 Mq and m2 = m bh = 10 Mq, giv- 
ing fc 2PN = at x = 0.141042, for r/ = 0.204082 (m ns = 
2.0 Mq, m bh = 5 Mq) fc 2 p N = at x = 0.139319, and for 
77 = 0.058594 (m ns = 1.0 Mq, m bh = 15 Mq) fc 2 p N = 
at x = 0.141884. Similarly, for specific values of x, k^pN 
becomes zero. As Eqn. (j4|) is quadratic in powers of x, for 
each value of 77, there exist two values of x at which A;3pn 
becomes zero. As examples, for i) = 0.058594, &3pn = 
at x = 0.254246,0.012740, for 77 = 0.107725, fc 3 PN = at 
x = 0.248754, 0.013051, and for 77 = 0.204082, fc 3 PN = 
at x = 0.236489, 0.0139767. As we have already mentioned 



We will discuss more about these mass ranges later. 



Figure 1. Variation of R with x for an NSBH binary with 77 
values 0.204082 (solid line), 0.107725 (dashed line) and 0.058594 
(dotted line) respectively, keeping x ^ 0.59. 



that x cannot be less than 0.59, both &2pn an d ^pn are al- 
ways non-vanishing and positive, so they add up in the total 
fcpN- That is why, it is worthwhile to study the significance 
of &3pn in comparison with &2pn. 

Expressions of fciPN and A;2pn in terms of basic or- 
bital elements like orbital periods, eccentricities and masses 
of the component s are simple and were first given by 
iDamour fc SchafeJ ij 19881 . hereafter DS88). On the other 
hand, expressing &3PN in terms of only those elements is 
rather difficult. As an alternative simple approach, we de- 
cide to compute the ratio 



R = 



fc3PN/fc2PK 



1 



fepN/fciPN ~ 24 [-5 + 277 + (35 - 10r/)x] 2 

[(120 - 12077 + 96?7 2 ) + 
(-10080 + 13952?7 - 1237r 2 77 - 1440t7 2 )x 
+ (36960 - 4000077 + 615tt 2 77 + 168077 2 )x 2 ] , 



(5) 



which depends only on 77 and x. In Fig[T]we plot R against x, 
x being in the range of 0.59 — 100, for the above mentioned 
values of 77. We see that R ~ 1.2, implying fcspr-j/fepN — 

fc2PN/fclPN. 

fciPN and &2pn are given in terms of orbital elements in 
DS88 as: 

3/3 2 



fciPN = 



1 



where 
and 

h = 



/3o 



3/?o 4 fo 



(GMn ) 1/3 



(6) 



(8) 



39 



27 



13 2 ,1 2 ,13 



(9) 



Here e t is t he proper time eccentricity (Da mour fc Deruellei 
1 19851 . 1 19861 ) which becomes the ordinary eccentricity in case 
of Newtonian dynamics, no = 2n/P or b, P or b is the orbital 
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Figure 2. Variation of t^ipH with orbital periods and orbital ec- 
centricities for an NSBH binary having m ns = 1.4 Mq, mth = 
10 Mq . X and Y axes represent orbital periods in hours and or- 
bital eccentricities respectively while the color code represents the 
values of i^ipn in cleg yr _1 . The contours of t^ipN = 100, 10, and 
1 deg yr _1 are also shown from left to right. 
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Figure 3. Variation of A^pn/^ipn with orbital periods and or- 
bital eccentricities for an NSBH binary having m n s = 1.4 Mq, 
m bh = 10 Mq . X and Y axes represent orbital periods in hours 
and orbital eccentricities respectively while the color code rep- 
resents the values of A^pn/^ipn- The con t° urs of ^pn/^ipn = 
10 -3 , 10 -4 , and 10~ 5 are also shown from left to right. 



period of the binary, M = m ns +m bh , x ns = m n3 /M, x bh = 
m bh /M, and ^i, 2 pn = n (fciPN + fera). 

For double neutron star (DNS) systems, observed ui 
is modelled using the 1PN term only, which is justified, 
as even for the most relativistic system PSR J0737-3039, 
fepN/fciPN = 2.6 x 1CT 5 , but the 2PN term should be 
considered if a more relativistic DNS system is discovered. 
The observed value of lj (16.899 deg yr _1 ) for PSR J0737- 
3039 l|Kramer et fflZ.l l2006) is the maximum among all known 
DNSs. For the sake of comparison, in Fig. [2] we show the 
variation of Wipn = no A;ipn with orbital periods and ec- 
centricities for an NSBH binary having m n3 = 1.4 Mq, 
rribh = 10 Mq. The X and Y axes represent orbital periods in 
hours and orbital eccentricities respectively while the color 
code represents the values of ojipn in deg yr _1 . The contours 
of wipn = 100, 10, and 1 deg yr _1 are also shown from left 
to right. Here one should remember that for a relativistic bi- 
nary, both P or b and et decrease with time due to the emission 
of gravitational waves. As an example, for an NSBH binary 
having m ns = 1.4 Mq, m bh — 10 Mq, P orb — 10 hours and 
e t = 0.3, 1PN values ar e P orb = -1 .04 x 10" 12 s s^ 1 and 
e t = -6.73 x 10 -18 s _1 (|Petersll 19641 ). 

In Fig [3] we plot fepN/fciPN f° r the same NSBH bi- 
nary. The X and Y axes have the same meaning and the 
color code represents the values of fepN/feiPN- The con- 
tours of A;2pn/A;ipn — 10 -3 , 10 -4 and 10 -5 are also shown 
from left to right. As the phase-space with very high ec- 
centricities and very short orbital periods is unlikeljQ we 
can say ^pn/^ipn < 10 -3 . As we have already seen 
that &3PN/&2PN ~ fepN/^iPN, it can be concluded that 
fc3PN/fciPN < 10~ 6 . Since the accuracy as good as 10~ 5 in 
the measurement of (li has been already achieved for many 
binary pulsars, one might aim for the same for NSBH bina- 



2 Such systems are unlikely to form and even if they form un- 
der some extra-ordinary situation, they would merge very soon 
making the probability of detection very small. 



ries. In that case, one must consider 2PN effects, but 3PN 
and higher order PN effects can be excluded unless further 
accuracy is intended. 

Now, we wish to study the effect of the leading order 
spin-orbit coupling. According to DS88, 



where 



and 



3/3g 



(g s ,nsfis,ns + gs, b hPs, b h) ■ 



2tt 



c 2 a b h 



(10) 



(11) 



(12) 



Gm bh 

Here I na is the moment of inertia of the neutron star, a b h = 
(2nl b h) I \cm b h P a , b h) is the spin parameter of the black hole 
(Ibh is the moment of inertia of the black hole and P 3]b h is its 
spin period), Ps,ns is the spin period of the neutron star. It 
is well known that a bh ^ Gm bh /c 2 , with a bh>max = Gm bh /c 2 
giving the "maximally rotating black hole" . From Eqns ((6]l 
and (I10II. we can write 



W = 



sPs,ns + gs, b hfis,bh) Po ■ 



(13) 



felPN 

Let us denote s ns as the unit spin vector of the neutron 
star, Sbh as the unit spin vector of the black hole, k as the 
unit orbital angular momentum vector, h ns as the unit vec- 
tor from the earth towards the neutron star perpendicular 
to the plane of the sky, hbh as the unit vector from the earth 
towards the black hole perpendicular to the plane of the sky. 
Let us also assume that i is the inclination angle of the or- 
bital plane with respect to the sky plane. Then, following 
DS88, 
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and 



6(1 - e 2)i/2 s in- 



[(3 sin i — 1) k + cosi h„ 



(14) 



Xbh (iXbh + 3x„s) r,_ . 2. 1N , 

ffsM = — r: 2vT75 ■ o- 3 sm i - 1) k - 

6(1 — ej) 1 / 2 sm 2 j L 



hbh] • s b h (15) 



will have the maximum value when s ns is parallel 



to the vector (3 sin i — 1) k+ cos i h ns , giving 



1 1/2 



3 + 



(4x ns + 3x bh ) 

6(1 - e 2 ) 1 / 2 ■ 



(16) 



For any other orientation of the vectors, g 3 .ns would be dif- 
ferent, but it would always be in the range of +g s ,ns, max to 
—g s ,ns, max- In case of s ns || k, one obtains 

X na (Ax na + 3x bh ) 



3(l-e 2 )V2 



(17) 



Similarly, g s ,bh will have the maximum value when Sbh 
is parallel to the vector (3 sin 2 i — 1) k + cos i hbh, giving 



gs.bh, max 



3 + 



1/2 



Xbh {^x bh + 3x ns ) 
6(1 -e 2 ) 1 / 2 ' 



(18) 



For any other orientation of the vectors, g s ,bh would be dif- 
ferent, but it would always be in the range of +g s ,bh, max to 
-g s ,bh, max- In case of s b h || k, one obtains 



gs,bh, 



Xbh (4:X b h + 3x„ a ) 



3(1 



s 2Xl/ 2 



(19) 



/s.ns, max 



From the above equations, it is clear that g s 
differs from g s>ns> \\ (or g Slb h, max from g s>bhy y) only by the 



factor i [3 + 



,1/2 



, which lies between 1.7 — 1.0 for i in 



2 1 sin' 3 

the range of 20° - 90°. 

To study the spin-orbit coupling effects for an NSBH 
binary, one needs to know the moment of inertia of the neu- 
tron star and the spin parameter of the black hole in addi- 
tion to their masses, orbital parameters and the spin period 
of the neutron star. The value of the moment of inertia of 
a neutron star of known mass depends upon the Equation 
of State of the matter, which is not very well constrained 
at present. We choose the values of the moment of inertia 
as 0.8 x 10 38 , 1.0 x 10 38 , and 1.4 x 10 38 kg m 2 respectively 
for neutron stars of masses 1.0, 1.4 and 2.0 M§ accord- 
ing to th e predictions of most of the standard Equations 
of State fArnett fc Bowers! Il977l; iBeiger fc Haensell |2002| ; 



lLattimer fc Schutj2005l ; lBagchill201oT ) 



First, we compare the spin-orbit coupling due to the 
neutron star with that due to the black hole by computing 
the ratio 



Su = 



9s , Ti 



II & 



gs,bh, \\Ps 



(20) 



This ratio is independent of P or b and et- In Fig. [4] we plot Sy 
against P St „ s for the same NSBH binary for different values 
of a b h/abh. m.ax- Here we use the limit a bh = 0.998 a,bh,max 
as given bv lThornd { 1974). It is clear that usually the effect 
due to the spin of the black hole dominates, i.e. Sn < 1, 
unless either the neutron star is rotating very fast or the 
black hole is rotating very slowly. As some examples, to get 
Sn > 1, a b h must be less than 0.01 a b h, ma x for P s , ns = 4 ms, 



0,1 



Figure 4. Variation of Sn with P s ,ns for an NSBH binary having 



1.4 M, 



0, ™bh 



10 Mq. The lines are marked with the 



corresponding values of a bh /a bh , 



a b h must be less than 0.005 a b h,max for P s , ns = 10 ms, and 
a b h must be less than 0.001 a b h,max for P s ,ns = 40 ms. 

Now, to answer the question whether the total spin- 
orbit coupling effect is comparable to 1PN effect, we com- 
pute the value of Wn where we replace g s ,ns by g snS) \\ an d 
g s ,bh by g s ,bh, || in Eqn. (jT3j) . We plot W\\ against orbital pe- 
riods and orbital eccentricities for the same NSBH binary in 
Fig. [5] As expected, W\\ increases as either the black hole or 
the neutron star or both become faster. W\\ increases with 
the increase of et and decreases with the increase of P orb pro- 
vided all other parameters are unchanged. It is clear that 
usually Wn > 10 -5 unless the black hole is rotating very 
slowly, e.g. for a bh = 0.01 a bh:max , the whole P orb - e t space 



has W« > W~ b for both P s 



0.1 s and P s . 



0.01 



s (even for P a , n s = 0.001 s which has not been shown in 
the figure). For a b h = 0.001 a b h,max, a significant portion 
of P orb — et space has W\\ < 10 -5 when P a , ns = 0.1 s, but 
for P s ,ns = 0.01 s only the region with high P orb and small 
et has Wn < 10~ 5 . These facts imply that fcgo term might 
have significant contribution to the observed u) for an NSBH 
binary depending on the spin and orbital parameters of the 
system. But here we have assumed that the spin vectors of 
both the black hole and the neutron star are aligned with 
the orbital angular momentum vector. In reality, they can 
be misaligned, specially if the NSBH system was formed 
through an exchange interaction. For some specific orienta- 
tions of the spin vectors, combined with appropriate values 
of other relevant parameters, kg§ can be very small and even 
zero. But as only a small window in the total phase-space 
can lead to this situation, we consider this situation to be 
unlikely, although not impossible. 

In general, W\\ is not so high that we should expect the 
non-leading order of the spin-orbit coupling to be effective. 
As an extreme example, let us consider the NSBH system 
to have P or b = 0.5 hours, e t = 0.5, a b h = 0.9 a b h,max, and 
P s ,ns = 1 ms, giving Wn = 7.16 x 10~ 3 . If we assume that 
LO-SO/1PN - NLO-SO/LO-SO, then NLO-SO/1PN be- 
comes ~5x 10 . In this case, the non-leading order would 
contribute to the observed value of u. But this set of param- 
eters is not very likely, most of the parameter space would 
probably lead to NLO-SO/1PN < 10~ 6 . Further exploration 
is needed before making any strong statement, which is be- 
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Figure 5. Variation of Wn with orbital periods and orbital eccentricities for an NSBH binary having m ns = 1.4 Mq, m bbb = 10 Mq. X 
and Y axes represent orbital periods in hours and orbital eccentricities respectively while the color code represents the values of W\\ . The 
contours for W\\ = 10~ 5 also have been plotted, whenever in the range. 



yond the scope of this paper due to the lack of the analytical 
expression for kso°- 

So far, we have displayed the results only for a sample 
NSBH binary with the neutron star of mass 1.4 Mq and the 
black hole of mass 10 Mq except for the variation of R with 
x. This choice of the mass o f the black hole may seem to 
be rather high as I6zel et al\ <|2010h found that the masses 
of the known black holes in low mass X-ray binaries have 
a narrow distribution at ~ 7.8 Mp, population synthesis 
for N SBH binaries by IPfahl. Podsiadlowski. fc Rappaportl 
also predicted the mean mass o f the black hole s 
to be around 7 Mq. On the other hand, iFarr et al\ (120111 ) 
found that the inclusion of high mass systems resulted the 
mean of the masses of known b lack holes to be 10-11 Mq. 
iFaucher-Giguere fc Loebl l|201ll ) used the black hole mass to 
be 10 Mq in their study of the exchange scenario. More- 
over, the higher value of the mass of the black hole makes 
the system more relativistic. Nevertheless, we have checked 
that our main conclusions remain unchanged for realistic 
ranges of the masses, i.e. the neutron star mass in the range 
of 1-2 Mq and the black hole mass in the range of 5-15 
Mq. We have already seen that, R is always around 1.2 for 
the above ranges of masses. The effects of the variation of 



the mass of the neutron star on fepjv/fciPiv, S\\ and Wu 
are insignificant, but the variation of the mass of the black 
hole has small but significant effects. As an example, for 
m ns ~ 1-4 Mq, P or b — 10 hours, e t = 0.3, the values of 
k2PN/k 1PN are 2.1 x 10" 5 , 3.2 x 10" 5 , and 4.1 x 10~ 5 for 
rribh = 5, 10, 15 Mq respectively. The fact that the phase- 
space for k,2PN /fcipjv > 10 -3 is unphysical, is always true, 
and thus our main conclusion that k3PN /feipjv < 10 -6 re- 
mains unaltered. Similarly, the maximum value of Su in Fig. 
[4] is 41 where P Si „ s = 0.001 s and a bh /a bh , ma x = 0.001 for 
rribh = 10 Mq, which becomes 87 if mbh = 5 Mq and 27 
if m bh = 15 Mq keeping P s , ns and a bh /a bht7nax fixed. The 
minimum value of Sy in Fig.[4]is 4.1 x 10~ 5 where P s , ns = 1 
s and a b h/a b h,max = 0.998 for m b h = 10 Mq, which becomes 
8.7 x 10" 5 if m bh = 5 Mq and 2.7 x 10~ 5 if m bh = 15 Mq 
keeping P 3 , n3 and a bh /a bh}max fixed. For m ns = 1.4 Mq, 
P orb = 10 hours, et = 0.3, P s ,„ s = 10 ms, a bh = 0.3 a bh , m ax 
the values of W\\ are 5.6 x 10~ 4 , 7.7 x 10" 4 , and 9.2 x 10~ 4 
for m b h = 5, 10, 15 Mq respectively. 

We have considered only the short period NSBH bina- 
ries, although very wide (P orb > 10 00 days) NSBH systems 
can also form (|Kiel fc Hurlevl 2009). We have excluded such 
wide systems from our study as it is clear that higher or- 
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der effects would be negligible for such systems. For com- 
pact NSBH systems, our chosen ranges for P or b and et are 
sufficiently wide to incorporate different formation scenar- 
ios, e.g. normal evolutionary formation and the formation 
through exchange interactions. The same logic applies for 
the choice of wide ranges for P s ,ns and ath- 



3 DISCUSSIONS 

In this study we have not dis cussed the e ffect of the spin- 
spin interaction. According to IWexl (|f 995h . spin-spin inter- 
action is significantly smaller than 2PN effect, and can be 
neglected. 

Five post-Keplerian parameters can be measured for a 
relativistic binary, the rate of the advance of the perias- 
tron (lu), the rate of change of the orbital period (P or b), 
Einstein delay parameter (7), Shapiro ra nge parameter (r), 
and Shapiro shape parameter (s = sini) l)Damour fc Tavlorl 
1 19921 ; II orimcr fc Kramer|[2005l i. We have discussed ill in de- 
tails, and briefly mentioned about IPN P or b- Other effects 
might also contribute significantly to the observed value of 
Porb, like the tail effect an d the s pin-orbit coupling as dis- 
cussed by iRieth fc Schafer] (]f997). Another interesting hy- 
pothesis is the outspiral of the NSBH binary due to the 
evaporation of the black hole caused by the e xtra spatial 
dimensions. According to lSimonetti et all (|201lT ) . this effect 
will be measurable for an NSBH binary and will test the 
validity of extra dimensions. Lowest order values of 7 and r 
for an NSBH system having m ns — 1.4 Mq, iribh ~ 10 Mq, 
Porb = 10 hours and e t = 0.3 are 0.013 s and 4.9 x 10~ 5 s. 
These are orders of magni tude higher than those for J0737- 
3039 (Kr amer et ffll.|[2006T ). So it might be worth exploring 
higher order effects for these parameters in case of NSBH 
binaries. 

In addition to contributing to the periastron advance, 
spin-orbit coupling also leads to the geodetic precession 
where both spin vectors precess around the total angular 
momentum of the system causing changes in the pulse pro- 
file. This effect h as been observed for some relativistic sys- 
tems {stairs 20 031) . The predicted rate o f precession upto the 
leading order IjLorimer fc Kramerll2005l ) for an NSBH binary 
having m ns — 1.4 Mq, nibh = 10 Mq, P or b = 10 hours, and 
et = 0.3 is 2.17 deg yr _1 , which is of the same order as 
expected from DNS systems, e.g. smaller than that of PSR 
J0737-3039 (4.78 deg yr _1 ), but larger than PSR B1913+16 
(1.21 deg yr" 1 ). This precession contributes to the observed 
change in the pr ojected semi-major axis (x = a p sin i/c) and 
the eccentricity ((Lorimer fc Kramer! 1200a ). This change in 
x is an addition to the possible change caused through the 
change of i caused directly by the spin-orbit coupling (Eqn. 
5.15ofDS88). 

Several other manifestations of general relativity 
are effective in case of NSBH b inaries, like frame- 
dragging effect (|Wex fc Kopei kin 1999) and tidal interaction 
|Daniour fc Nagai 20ld . and references therein), etc. Study 
of these effects was not the aim of this paper. 



4 CONCLUSIONS 

In this work, we have shown that the 2PN and the leading 
order spin-orbit coupling terms can have significant contri- 
butions to the observed value of the periastron advance of an 
NSBH binary, while the 3PN and higher order terms can be 
neglected. It will not be difficult to incorporate &2PN term in 
the future timing code, as this has a relatively simpler form. 
On the other hand, kg§ might be difficult to incorporate as 
it involves few difficult-to-know terms like the moment of 
inertia of the neutron star, spin of the black hole and orien- 
tations of the vectors. Moreover, higher order terms will be 
important for modelling gravitational wave emission from 
such systems. So further theoretical study in this field is 
necessary before the discovery of the first NSBH binary. 
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